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Abstract 



Consider the general scalar balance law dtu + Div/(/;, x, u) = F{t, x, u) in several space 
dimensions. The aim of this note is to estimate the dependence of its solutions from 
the flow / and from the source F . To this aim, a bound on the total variation in the 
space variables of the solution is obtained. This result is then applied to obtain well 
posedness and stability estimates for a balance law with a non local source. 
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1 Introduction 

The Cauchy problem for a scalar balance law in N space dimension 



is well known to admit a unique weak entropy solution, as proved in the classical result 
by Kruzkov [12', Theorem 5]. The same paper also provides the basic stability estimate on 
the dependence of solutions from the initial data, see [12i Theorem 1]. In the same setting 
established in [1^, we provide here an estimate on the dependence of the solutions to (|1.1|) 
from the flow /, from the source F and recover the known estimate on the dependence 
from the initial datum Uq- A key intermediate result is a bound on the total variation of 
the solution to (II. ip . which we provide in Theorem 12.51 

In the case of a conservation law, i.e. -F = 0, and with a flow / independent from i, x, 
the dependence of the solution from / was already considered in [3], where also other results 
were presented. In this case, the TV bound is obvious, since TV {u{t)) < TV(no). The 
estimate provided by Theorem l2.5l slightly improves the analogous result in [3", Theorem 3.1] 
(that was already known, see OdB]), which reads (for a suitable absolute constant C) 




(1.1) 



u{t) - v{t) 



Li(R'V;]R) 



< \\Uo - Vo 



Li(RiV;iR) +CTV(mo) Lip(/-g)t. 
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Our result, given by Theorem l2.6t reduces to this inequality when / and g are not dependent 
on t, X and F = G = 0, but with C = 1. 

An flow dependent also on x was considered in [11[9], though in the special case f{x, u) = 
l{x) g{u), but with a source term containing a possibly degenerate parabolic operator. 
There, estimates on the L"*" distance between solutions in terms of the distance between 
the flows were obtained, but dependent from an a priori unknown bound on TV . 
Here, with no parabolic operators in the source term, we provide fully explicit bounds both 
on TV {u{t)) and on the distance between solutions. Indeed, remark that with no specific 
assumptions on the flow, TV may well blow up to +00 at t = 0+, as in the simple 

case f{x,u) = cosx with zero initial datum. 

Both the total variation and the stability estimates proved below turn out to be optimal 
in some simple cases, in which optimal estimates are known. 

As an example of a possible application, we consider in Section [3] a toy model for a 
radiating gas. This system was already considered in [5l [HI [ini [HI IS [HI [151 [T7] . It consists 
of a balance law of the type (jl.ip . but with a source that contains also a non local term, due 
to the convolution of the unknown with a suitable kernel. Thanks to the present results, 
we prove the well posedness of the model extending [HI Theorem 2.4] to more general flows, 
sources and convolution kernels. Stability and total variation estimates are also provided. 

This paper is organized as follows: in Section [21 we introduce the notation, state the 
main results and compare them with those found in the literature. Section [3] is devoted to 
an application to a radiating gas model. Finally, in sections [H and [5] the detailed proofs of 
theorems 12.51 and 12.61 are provided. 



2 Notation and Main Results 

Denote M+ = [0, +oo[ and M+ = ]0, +oo[. Below, iV is a positive integer, O = M+ x x M, 
B{x,r) denotes the ball in with center x G and radius r > 0. The volume of the 
unit ball 5(0,1) is ujn. For notational simplicity, we set ujq = 1. The following relation 
can be proved using the expression of a; at in terms of the Wallis integral Wn- 



2Wn where Wn = / (cos0)^d0. (2.1) 

Jo 



In the present work, 1^ is the characteristic function of the set A and 6t is the Dirac 
measure centered at t. Besides, for a vector valued function / = f{x,u) with u = u{x), 
Div/ stands for the total divergence. On the other hand, div/, respectively V/, denotes 
the partial divergence, respectively gradient, with respect to the space variables. Moreover, 
du and dt are the usual partial derivatives. Thus, Div/ = div/ + duf • Vn. 

Recall the definition of weak entropy solution to (jl.ip . see |12^ Definition 1]. 

Definition 2.1 A function u E L°°(M+ x M^;]R) is a weak entropy solution to lll.l]) if: 
1. for any constant A; € M and any test function ip G C^(M+ x M^;M_|_) 



{u - k) dt(p + (/(t, X, u) - f{t, X, k)) ■ Vip + {F{t, x, u) - div f{t, x, k)) cp 
xsign(u — k) dx dt > 0; 



(2.2) 



2. there exists a set E of zero measure in M_f_ such that for t G 
is defined almost everywhere in and for any r > 



\ £ the function u{t, x) 



lim 

t-*o,tm+\e JBio,r) 



\u{t, x) — Uo{x)\dx = . 



(2.3) 



Throughout this paper, we refer to [HUB] as general references for the theory of BV func- 
tions. In particular, recall the following basic definition, see [U Definition 3.4 and Theo- 
rem 3.6]. 



Definition 2.2 Let u G LJ^^ 



;M). Define 
TV(u) = sup < / udivipdx: 

u G Lii„^(M^;M) : TY{u) < +cx)} . 



and 



< 1 



BV( 



The following sets of assumptions will be of use below. 



(HI) 



(H2) <^ 



(H3) <^ 



/ G C^{Q;R^) 
c>„/G L-(0;M^) 
a„(F-div/) gL~(J^; 
/ G C^{Q;R^) 

Vduf G L~(17;M^x^) 

aA/GL~(17;IR^) 

dtdivf G L~(J7;M) 
/ G C^{n;R^) 

a„/GL-(f];M^) 



F G C^(17;IR) 

I) F-div/ G L°°(1^;M) 
F G Ci(17;M) 

/ / ||V(F-div/)(t,x,.)|l 
dtF G L°°(0;R) 



L°°( 



dx dt < +00 



F G C°(J1; M) duF G L°°(f]; R) 
\\{F -divf){t,x,-)\\^^, . 



dx dt < +00 



The quantity F — divf has a particular role, since it behaves as the "true" source, see (j2.6p . 
We note there that the assumptions above can be significantly softened in various specific 
situations. For instance, the requirement that / be Lipschitz, which is however a standard 
hypothesis, see [21 Paragraph 3], can be relaxed to / locally Lipschitz in the case / = 
f{u) and F = 0, thanks to the maximum principle [121 Theorem 3]. Furthermore, the 
assumptions above can be obviously weakened when aiming at estimates on bounded time 
intervals. 

Assumptions (HI) are those used in the classical results |12^ Theorem 1 and Theo- 
rem 5]. However, we stress that the proofs below need less regularity. As in [12], we 
remark that no derivative of / or F in time is ever needed. Furthermore, / needs not be 
twice differentiable in u, for the only second derivatives required are Vxduf and V^/. 

We recall below the classical result by Kruzkov. 



Theorem 2.3 (Kruzkov) Let (HI) hold. Then, for any Uo G L°°(M^;M), there exists 
a unique weak entropy solution u to in L°° ^M+; Lj^^(R^; M)^ continuous from the 

right. Moreover, if a sequence G L°°(M^;M) converges to Uq in LJ^^^, then for all t > 
the corresponding solutions u^{t) converge to u(t) in Lj-^^. 



Remark 2.4 Under the conditions (H2) and Jjg^/[gjv||(-^ — div/)(t, x, dxdt < 

+00, see (H3), the estimate provided by Theorem 12.51 below, ahows to use the technique 
described in [71 Theorem 4.3.1], proving the continuity in time of the solution, so that 

ueC^ (m+;L1,JM^;M)). 

2.1 Estimate on the Total Variation 

Recall that [9, Theorem 1.3] and 01 Theorem 3.2] provide stability bounds on p.ip . in the 
more general case with a degenerate parabolic source, but assuming a priori bounds on 
the total variation of solutions. Our first result provides these bounds. 

Theorem 2.5 Assume that (HI) and (H2) hold. Let Uo G BV(M^;M). Then, the weak 
entropy solution u of satisfies u{t) E BV(R^;M) for all t > 0. Moreover, let 

Ko = NWn ((27V + 1) ||Va„/||L.o(f,.iRivx^) + \\duF\\^^^^.^^ (2.4) 
with Wn as in 112.1]) . Then, for all T > 0, 

TY {u{T)) <TV{uo)e^°^ + NWn [ e''°(^-*)/ \\V{F-divf){t,x,-)\\roodxdt. (2.5) 

Jo Jrn 

This estimate is optimal in the following situations: 

1. If / is independent from x and F = 0, then Kq = and the integrand in the right 
hand side above vanishes. Hence, (|2.5|) reduces to the well known optimal bound 
TV (u(t)) < TV(no). 

2. In the ID case, if / and F are both independent from t and u, then Kq = and (II. 1|) 
reduces to the ordinary differential equation dtu = F — div/. In this case, (|2.5|) 
becomes 

TV {u{t)) < TV(uo) + t TV(F - div/) . (2.6) 

3. If / = and F = F{t) then, trivially, TV (u(t)) = TV(uo) and ([23]) is optimal. 
A simpler but slightly weaker form of (j2.5p is 

TV (n(T)) < TV(no) e^°^ + NWn ^ ° ~ sup / \\V{F - dwf){t, x, -Ml dx 

te[o,T] Jm 

when the right hand side is bounded. 



2.2 Stabihty of Solutions with Respect to Flow and Source 

Consider now (jl.ip together with the analogous problem 

dtv + Bwg(t,x,v) =G{t,x,v) {t,x) G M+ x 
v{0,x) =Vo{x) X e M^. ^ ■ -* 

We aim at estimates for the difference u — v between the solutions in terms of f — g, 
F — G and Uq — Vq- Estimates of this type were derived by Bouchut &: Perthame in [3] 
when /, g depend only on u and F = G = 0. Here, we generalize their result adding the 
(i, x)-dependence. The present technique is essentially based on Theorem 12.51 



Theorem 2.6 Let {f,F), {g,G) verify (HI), {f,F) verify (H2) and {f -g,F-G) ver 
ify (H3). Let Uo^Vo^ BV(M^;M). We denote Ko as in 1^2.4^ and introduce 

K = 2iV||V(9„/||Loo(Q.iRiVxiV) + ||<9ni^|lL°o{Q;R) + ||t^n(-^ - 1 1 ^oo (q-jr) and M = ||<9m5|Il=°{Q; 

Then, for any T,R > and Xq G M^, the following estimate holds: 



+ 



x—Xo\\^R 



|ti(r, x) — f (r, x) Idx < e 



\uo{x) — Vo{x) I dx 



TVK)||a„(/-5)| 



+ NWn 



T ^Ko(T-t) _ ^K(T-t) 



|V(F-div/)(t,x, 



,dxdt \\duif - g)\ 







x-Xo\\<R+M(T-t) 



((i7-G)-div(/-5)) (t,x,-) 



dx dt . 



The above inequahty is undefined for k = Kq and, in this case, it reduces to (j5.17p . This 
bound is optimal in the fohowing situations, where Uo,Vo € L-'-(R^;M). 

1. In the standard case of a conservation law, i.e. when F = G = and f,g are 
independent of x, we have = and the result of Theorem ED becomes, see O 

Theorem 2.1], 



|u(T) ^'(^)||Li(]RiV;R) 



< \\Un 



+ TTY{uo)\\du{f -g)\ 



L°°(f7; 



2. If duf = dug = and duF = duG = 0, then Kq = k = Q and Theorem 12.61 now reads 



|n(r)-^;(r)|L, 



< \\uo - Vo\ 



+ 



[(i7_G)-div(/-<7)] (t) 



Li(M^;M) 

3. If {f,F) and {g,G) are dependent only on x, then Theorem 12.61 reduces to 



dt. 



\uiT)-viT)\L., 



< \\Uo - Vol 



+ T\\iF-G)-divif-g)\U 



The estimate obtained in Theorem 12.61 shows also that, depending on the properties of 
specific applications, the regularity requirement / G C^(r2; M'^) can be significantly relaxed. 
For instance, in the case f{t,x,u) = q{u) v{x) considered in [UIH], asking q of class C-*- and 
V of class is sufficient. See also Section [3] for a case in which the required regularity in 
time can be reduced. 

In the case of conservations laws, i.e. when F = G = 0, one proves that k < Kq and the 
estimate in Theorem 12.61 takes the somewhat simpler form 



/ |n(r,x) -?;(T,x)|dx < e''^ / 

J\\x — x„\\<R J \\: 



+ re''°^TVK)||a„(/-<7)| 



x-Xo\\<R+MT 



\uo{x) — Vo{x) \ dx 



Loo 



NWnT^c^"^ sup ( [ ||Vdiv/(t,x,-)||Loodx) \\du{f-g)\\ 



+ 

+ Te'^"'^ 



Loo 



sup 



|div(/-5r)(t,x, •)! 



dx 



te[0,T] J\\x-Xo\\<R+M(T-t) 

when the right hand side is bounded. In the case considered in [3l Theorem 3.1], / = f{u) 
Ko = and we obtain [3l formula (3.2)] with 1 instead of the constant C therein. 



u 



3 Application to a Radiating Gas Model 

The following balance law is a toy model inspired by Euler equations for radiating gases: 

dtu + Divf{t,x,u) = —u + K *x u . (3.1) 

It has been extensively studied in the literature when / = f{u), see for instance |10^ [TT| 
[131 [13 [H! for the scalar ID case, O [E] for ID systems, |8j for the scalar A^D case. 

The estimate provided by Theorem 12.61 allows us to present an alternative proof of 
the well posedness of (|3.ip proved in [8]. Furthermore, we add stability estimates on the 
dependence of the solution from / and K, in the case of / dependent also on t, x and with 
more general source terms. 

Theorem 3.1 Let {f,F) satisfy (HI), (H2) and (H3). Assume that 

t2r^TOo\/Tni -.^ TuiN , Tni\ ; t^' ^ t co f ttsi . -\-\t2 ,1 /m N , tt 



(K) K € (C^ n L°°)(M+ X ; M) and K G L°° (^M+; W 
Then, for any Uq € (BV H L-'-)(M^;M), the Cauchy problem 

I dtu + Divf{t,x,u) = F{t,x,u) + K {t,x) G M+xM^ 



u(0, x) = Uo{x) X G 



TV (3.2) 



admits a unique weak entropy solution u G C° ^]R+; L-'-(M^; M) j . Moreover, denoting 
k = ||-^|Il°o{r_,. li(r^ r))' f'^''" T > 0, the following estimate holds: 

TV(u(r)) < 6^'^°+^^^'=)^ TV(uo) 

+iVTy7V /"^ e('^°+^^^'=)(^-*) [ ||V(F-div/)(t,x,-)||Toodxdt. 
Jo Jrn 

If F{t, X, 0) - div/(t, X, 0) = for all t G [0, T] and x G M^, then 

2. Let K satisfy (K) and call u the solution to k3.2\) with K replaced by K . Then, 

\\u{T)-u{T)\\ <\\uo\\^,^^^N.^)—-^ K-K (3.3) 



k-k 



Proof. Fix a positive T (to be specified below) and consider the Banach space X = 
CO ([0,r];Li(M^;M)) equipped with the usual norm \\u\\yr = WuW^^aom ■l^CR^-r')')- Define 
on X the map T so that T{w) = n if and only if u solves 

{dtu + Ti\Yf{t,x^u) = F{t^x^u) + K*xW {t,x) G 
ti(0, x) = Uo{x) X G 



in the sense of Definition 12.11 Note that the source term does not have the regularity 
required in (HI). However, by the estimate in Theorem 12. 6^ we can prove that ()3.4p does 
indeed have a unique weak entropy solution, see Lemma 13.21 for the details. The fixed 



I 



points of T are the solutions to . By Theorem 12.31 and Remark \2.4:\ Tw € X for all 
w & X. We now show that T is a contraction, provided T is sufficiently small. Note that 

Ko = NWn{{2N + l)\\Vduf\\^^ + \\duF\\^^) 

K = 2N\\Vduf\\^^ + \\duF\\^^. 

Moreover, by Theorem 12.61 



d{Twi,Tw2) = sup \\Twi-Tw2\\i^i 

te[o,T] 

( e'^* - I \ 

< sup sup \\K{t) *^ {wi - W2){t)\\^j_ 

te[o,T] \ ^ Te[o,i] J 

< sup ||K(r)||j^i - u;2)(t)||li 

^ rG[0,T] 

e^^ - 1 

< kd{wi,W2)- 

K 

Therefore, T is a contraction as soon as T is smaller than a threshold that depends only on 
ll^«-^llL-(n:R)' l|V9„/||Loo(n:MivxiV) and on ||i^||L-(R+;Li(R^;iR))- Therefore, we proved the 



well posedness of (j3.2p globally in time. 

Consider the bound on TV By Theorem 1 2. 5 1 

TV {u{T)) < TV(uo) + NWn £ e^"^^-'^ j^^ ||V(F - div/)(t,a;, •)||Loo(M;MiV) dxdt 

+NWn [ e'^°(^-*)A;TV(n(t)) dt 
Jo 

and an application of Gronwall Lemma gives the desired bound. 

We estimate the L-*- norm of the solution to (j3.2p . comparing it with the solution to 

j dtu + Bwf{t,x,u) = F{t,x,u) + K {t,x) G M+ x 

I n(0,x) = X G R^. ^ ^ 

By assumption, solves ()3.5p . hence it is its unique solution. Then, evaluating the distance 
between the solutions of (13. 2|) and (13. 5p by means of Theorem 12.61 we get 

e~'^'^\\u{T)\\^^^^^,^^ < ||no||Li(]Riv.B;) + ^ e'""^ j^^\K *^u{t,x)\dxdt 
and, thanks to Gronwall Lemma, we obtain: 



The final estimate p.3p follows from Theorem 12.6 



e '^'^Wiu-u 



< 


K 


-k 


< 


K 


-k 









rT 

L-(M+;Li(RiV.M)) Jo ^~''*II"(*)IIl1(R'V;R) + e-'^'Wiu - u) (t) (rN dt 

and thanks to Gronwall Lemma, we get the result. 

The continuity in time is proved as described in Remark 12.41 □ 



Lemma 3.2 Let f,F satisfy (HI) and K satisfy (K). If w £ L°°{R+ x M^;M), then the 
estimates in Theorem \2.5\ and in Theorem \2.(A apply also to \3.4^. 



Proof. Fix positive T, R and let Wn be a sequence of C°° functions converging to w in 
L-*- ( [0, T] X ]R^;M). Apply Theorem 12.31 to the approximate problem 



\ dtu + V)iYf{t,x,u) = F{t,x,u)+K*^Wn it,x) G R+xM^ 

I u{0,x) = uo{x) X e ^ ^ 

to ensure the existence of its weak entropy solution n„. Apply Theorem 12.61 to estimate 
the distance between Un and ii„_i: 

- M„_i||loc([o_T];L1(M^;E)) ^ / e''^^"*) / \K * {Wn - Wn^l){t,x)\dxdt 





< e^'^ k \\Wn — Wn- 



1|Ili([0,T]xI[ 



showing that the form a Cauchy sequence. Their limit u solves (13. 2p . as it follows 
passing to the limit over n in the integral conditions (|2.2p - (j2.3p and applying the Dominated 
Convergence Theorem. The estimates in theorems 12.51 and 12.61 are extended similarly. □ 



4 Proof of Theorem 12.5 



Lemma 4.1 Fix a function //i G C^(M_|_;M+) with 

supp(;ui) C [0,1[, / r^-Vi(Odr = -^— , fi[<0, = for n > I. (4.1) 
Define 

M^) = wwf^V (4.2) 



Then, recalling that ujq = 1, 



H{x)dx = 1, (4.3) 
(||x||) dx = ^^^^^ / ||x||^i (||x||) dx, (4.4) 
/ |lx|l ||V^(x)|| dx = - |lx||/i'i (||x||) dx = N, (4.5) 



[ \\xf fi[ {\\x\\) dx = -(iV + 1) / ||x||^i (||x|l) dx. 



(4.6) 



Proof. The first relation is immediate. Equalities (14. 5p and (14. 6p follow directly from an 
integration by parts. Consider ()4.4|) . The cases = 1, 2, 3 follow from direct computations. 
Let A^ > 4 and pass to spherical coordinates {p,6i, . . . , 6n-i), 



xi = p cos 9n-i 



X2 = p sm9]\f^i cos6n-2 



xn-i = p smt)N-i smt>]\[-2 ■ ■ ■ costii 
xj\f = p sin0Ar_i sm07v-2 • • • sin^i 



with p G M+, 6*1 G [0, 27r[ and Oj € [0, vr] for j = 2, . . . , - 1. If iV > 4 

jxil pi (||x||) dx 



f r ... r icoseN-ilp"" pi{p) iT\{sm9jy-A deN-ideN-2---d9idp 

Jk+ Jo Jo Jo \ I 

/ [cos6'iv_i| (sin6'Ar_i)^"^ d6'Ar_i ) / p^pi{p)dp 
Jo / Jr+ 



= (jV-lViv-i .,^ ^ J I ||x||/ii(||x||) dx 
iV - 1 Nlvn Jrn 

= — ^^^^ / ||x||^i (||x||) dx 

completing the proof. □ 

Recah the foUowing theorem (see [TJ Theorem 3.9 and Remark 3.10]): 

Theorem 4.2 Let u € L/^^(M^;M), then u € BV(]R^;R) if and only if there exists a 
sequence Un in C°°(]R^;M) converging to u in LJ^^ and satisfying 

hm / II Vn„(x) II dx = L with L < oo . 

n^+oo J^N 

Moreover, TV(u) is the least constant L for which there exists a sequence as above. 

Proposition 4.3 Fix pi as in Let u € LJ-q^(M^;M) admit a constant C such that 

for all positive A, R and with p as in ^.S^ 

— I I |n(x) — n(x — I dx dz < C. (4-7) 

A Jrn Jb{xo,r) 

Then, u G BV(M^;R) and TV{u) < C/Ci, where 

Ci= [ |xi|/zi(||x|l) dx. (4.8) 



Note that Ci £ ]0, 1[. If moreover u G Ci(M^;M), then 

TV(n) = — hm — / / |n(x) — n(x — z)| ^(z) dxdz . (4-9) 



lU 



Proof. We introduce now a regularisation ofu: Uh = u*iJ,h, with jihix) = fJ-i (||a;||//i) /h^ . 
Note that Uh € C°°(M^;M) and Uh converges to u in LJ-^^, as /i — >■ 0. Furthermore, for R 
and /i positive, we have 



1 
A 
1 

< - 
- A 



B(xo,R) 
B{xo,R+h) 



Uh{x) — Uh{x — z) \ fi{z) dxdz 
u{x) — u{x — z) I ^{z) dx dz 



and 



< C 

Uhjx) -Uhjx - Xz) 
A 



Vuh{x — Xsz) ■ zds . 



Thanks to the Dominated Convergence Theorem, at the hmit A — s- we get 

I Vu/i(x) • z] /iidl^ll) dx dz < C . 

B{xo,R) 

Remark that for fixed x € B{xo, R), when Vuh{x) ^ 0, the scalar product Vuh{x) ■ z is 
positive (respectively, negative) when z is in a half-space, say (respectively, H~). We 



can write z = a- 



\\7uh{x) 

\Vuh{x) ■ z\ fJ,i{\\z\\) dz 



+ w, with a G M and w in the hyperplane H° = Vuh{x) . Hence 



Vuh{x) ■ z fii{\\z\\)dz + I Vuh{x) ■ {-z) fii{\\z\\)dz 
Vuh{x) ■ ziJLi{\\z\\) dz 



H° 



a ||Vu/i(x)|| iJii{\J + lltfll^) dwda 



\a\ II Vn/i(x)|| p,i{\J + \\w\\^) dw da 

= \\Vuh{x)\\ / |zi[ ;Ui(||z||) dz . 

Define Ci as in (j4.8p and note that Ci G ]0, 1[. Then we obtain, for all > 0, 

C 



Vuh{x)\\ dx < 

B(xo,R) 

Finally when i? — > oo we get f^^ ||VM/i(x)||drE < C/Ci and in the limit h — > 
Theorem 14.21 also TV(n) < C/Ci, concluding the proof of the first statement. 

Assume now that u G Ci(]R^;M). Then, using the same computations as above. 

1 



(4.10) 
0, by 



lim , 
A 



lim 



Ci TV(u) 



\u{x) — u{x — z) fi(z) dx dz 



Vu{x — Xsz) ■ z ds 



^i(||z||) dxdz 



completing the proof. 



□ 
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In the following proof, this property of any function u € BV(]R , M) will be of use: 

i{x)-u{x-z)\dx<\\z\\TY{u) forallzeM^. (4.11) 



), the general case will be 



For a proof, see [H Remark 3.25]. 

Proof of Theorem \2.5l Assume first that Uo S C^( 
considered only at the end of this proof. 

Let u be the weak entropy solution to (jl.ip . Denote u = u{t,x) and v = u{s,y) for 
{t, x), [s, y) G M+ X M^. Then, for all A;, / E M and for all test functions = ip{t, x, s, y) in 



Cl (M+ X 



,N\2.-- 



we have 



{u - k) dt(p + {f{t, X, u) - f{t, X, k)) Vxf + {F{t, x, u) - div f{t, x, k)) ip 

xsign(u — k)dxdt > 

for all {s,y) G M+ x M^, and 

{v - I) dsip + (/(s, y, v) - f{s, y, I)) Vyip + y, v) - div/(s, y, 

X sign(ti — /) dy ds > 



(4.12) 



(4.13) 



for all (t, x) G 



pN 



. Let G 



X M^;M+; 



, ^gC^ 



t>N.] 



) and set 



ip{t, X, s, y) = $(t, x) -^{t - s,x-y). 



(4.14) 



Observe that dtip + dg^p = ^dt^, Vx'P = ^ V^^> + ^> V^^*, Vytp = -^> V^.^'. Choose 
k = v{s,y) in (I4.12p and integrate with respect to {s,y). Analogously, take I = u{t,x) 
in (j4.13p and integrate with respect to (i,x). Summing the obtained inequalities, we get 



signm — V 



[u — V 



)^dt'^+ {fit, X, u) - fit, X, v)) ■ (V$) ^ 



+ if is, y, v) - f{s, y, u) - fit, X, v) + fit, X, u)) ■ (V^) cD (4.15) 



+ [Fit, X, u) — F(s, y, v) + div/(s, y, u) — div fit, x, v)) ip 



dxdt dyds > 0. 



Introduce a family of functions {Y^}^^q such that for any '& > 0: 



1 



Y^t) 



Y^is)ds 



V 



Y' G 

supp(y') c ]o,i 

y > 

Y'(s)ds = 1. 



— oo 

?? It? 



(4.16) 



Let M = ||5„/||loo(q.|5jv) and define for e,9,To,R > 0, Xq G M^, (see Figure [I]): 
xit) = Y,it)-Yeit-T) and ?/;(t, x) = 1 -Ye (||x - Xo|| - - M(ro - t)) > 0, (4.17) 



X 

1 



t 







1 






— > 



Oe T T + e Xo a b 

Figure 1: Graphs of x, left, and of V', right. Here a = R+M{To-t) and b = R+M{To-t)+e. 

where we also need the compatibility conditions To >T and Me < R+M{To — T). Observe 
that X ~^ l[o,T] and x' ^ ~ as e tends to 0. On x and we use the bounds 

X < l[0,T+£] and 'i-B{xo,R+M{To-t)) < fp < ^B(xo,R+M{To-t)+B) ■ 

In (I4.15p . choose = x{^) i'i^i^)- With this choice, we have 



dt^ = X' - M xYq and = -xYg 



Setting B{t, x, u, v) = \u — v|M + sign(u — v) [f{t, x, u) — f{t, x, v)) 
in ()4.15p becomes 



(4.18) 
the first line 



< 



{u-v)^ dt^ + {f{t, X, u) - f{t, x, v)) (V$) ^ sign(u - v)dx dt dy ds 

— v\x' — B{t, X, u, v)x Ye) dx dt dyds 

/ / / / \u — v\x d-xdtdyds 
Jr+ Jr'^ Jr+ Jr'^ 



since B{t,x,u,v) is positive for all {t,x,u,v) G 17 x M. Thanks to the above estimate and 
to (I4.15|) . we have 



> 0. 



{u — v)x' 

+ {f{s,y,v)-f{s,y,u)-f{t,x,v) + f(t,x,u))-{V^)<^ 
+ {F{t, X, u) - F(s, y, v) - div/(t, x, v) + div/(s, y, u)) cp 
xsign(u — v) dx dt dy ds 

Now, we aim at bounds for each term of this sum. Introduce the following notations: 

/ = / / / / \u — v\x' dxdtdyds , 
Jr^ Jr^ Jr^ 

{fit, y, v) - fit, y, u) + fit, X, u) - fit, x, v)) (VM-) $ 
X sign(ii — v) dx dt dy ds , 



Jt 





(/(s, y, v) - f{s, y, u) + f{t, y, u) - f{t, y, v)) {V^>) $ 
X sign(u — v) dx dt dy ds , 

[F{t, X, u) - F{t, y, v) - div/(t, x, v) + div/(t, y, u)) ip 
X sign(u — v) dx dt dy ds , 

{F{t, y, v) - F{s, y, v) - div/(t, y, u) + div/(s, y, u)) ip 
X sign(u — v) dx dt dy ds . 



Then, the above inequahty is rewritten as I + Jx + Jt + Lx + Lt > 0. Choose ^{t, x) = 
u{t) n{x) where, for r/, A > 0, /x e C^(M+; M+) satisfies and 

i/i(s)ds = l, z^i G C^(R;M+), supp(z^i) C 0[ . (4.19) 



We have 

/ < h + h 



where 



h = [ [ [ [ \u{t,x)-u{t,y)\ {Y^{t)-Y;{t-T)) ip^dxdtdyds, 
Jr+ Jrn Jr+ Jrn 

h = I I I I \u{t,y)-u{s,y)\{Yl{t)+Yl{t-T)) i,^dxdtdyds 

JW+ J'B.f^ JW+ JRJV 



and we get 



lim sup Ii < 

e~*0 



x-Xo\\<R+MTo- 



|n(0, x) — n(0, y)| iJ.{x — y) dxdy 



\u{T,x) - u{T,y)\ fi{x - y)dxdy , 



x-Xo\\<R+M{To-T) 

lim sup /2 < 2 sup / y) — n(s, y)| dy . 

e^O te{0,T}, J \\y-x^\\<R+\+M{To-t)+9 

sG]t,t+i7[ 

For Jx, we have that by (HI), / G C^(il;IR^) and therefore 

\\f{t,yiv) - f{t,y,u) + f{t,x,u) - f{t,x,v)\\ = 



V fl 




u JO 



Vduf [t, x{l — r) + ry, w) ■ {y — x) dr dw 



< l|V«9n/||Loo(Q.KivxiV)|Ix - y|| \u{s,y) - u{t,x)\ 



Then, using (j4.5|) 

Jx < !|V9„/||lo 

< llva„/|iLo 




||x — y|| \u{t, x) — u{s, y)| || V^*!! x dx dt dy di 

V 

\\x - y|| \u{t,y) - u{s,y)\ + \u{t,x) - u{t,y)\ 
V L 

x||V^'|| xV'dx dt dyds 



< N\\Vd.Jh^{T + e) sup / 

te[0,T+e], J\\y-xJ\ 



\u{t,y) - u{s,y)\dy 



se]t.t+vl 



+ l|Va„/||L. 




Jt < 



\\<R+X+M{Ta-t)^ 

X - y\\ \u(t,x) - u{t,y)\ \\Vi^\\dxdydt, 

IIV^'II $dx dt dyds 
< vWdtdufhoo [ [ [ [ \u{t,x)-uis,y)\\\V^\\<l>dxdtdyds. 

JR+ JR^ JK+ JR^ 




T+e 



JR^ J B{xo,R+M{To-t)+e) 

t fU 

dtduf{T,y,w) dw dr 




For Lx, we get 



Li 



Li + L2 where 



{dudivfit, x, w) + duF{t, y, w)) dw 



if sign(tt — v) dx dt dy ds, 



V{F — div/) (t, rx + (1 — r)y, u) ■ {x — y) dr 



xsign(u — v) dx dtdy ds. 
Then, recalhng (14.140 . the definitions ^ = ly fi, ^ = xi^, (fO) . (14.190 and (I4.17p . we obtain 



(T + e) sup / 

ie[0,T+s], J||y-Xo 



- ^^(•s,y)| dy 



+ 



se]t,t+ril 

T+e 



y-Xo\\<R+\+M(To~-t)+e 

u{t, x) — u{t, y) I — y) dx dy dt 



JR^ J\\x-Xo\\<R+]\I{To-t)+e 



L2 < V(F — div/) (t, y + r(x — y), u) \\x — y\\xi^ fJ-i^drdxdtdyds 

Jr+Jr'^Jr+ Jr'^Jo 

^ i f f ||V(F-div/)(t,y,-)||Lood2/dt) / ||x|l^(x)dx 
\Jo Jr^ J JR^ 

= AAfi [ ||V(F-div/)(t,y,.)||Lood?/dt 

Jo JRi^ 



where 

Ml = / ||x||^i (||x|[) dx. 

Concerning the latter term Lt 

Lt < v^N{R + MTo f {T + e){\\dtdiYf\\^^ + \\dtF\\j^^) 
Letting e, t/, ^ — > we get 



(4.20) 



lim sup Ii 




x-Xo\\<R+MTo 



\u{0,x) — u{0,y) \ fi{x — y) dxdy 




\x-Xo\\<R+M{To-T) 



I n(r, x) -u{T,y)\ii{x - y) dx dy , 



lim sup I2 = , 



limsupJi. < i|V(9„/||Loo / / 

£,-0,9^0 Jo JR' 



\x - y\\ \u{t,x) - u{t,y)\ 



lB(xo,R+M{To-t)) 
x||V/i(x — y)|| dxdy dt , 

lim sup Jt = , 

limsupLi < duf Wj^oo + \\duF\\j^oa) 

X / \u{t,x) -u{t,y)\fi{x-y)dxdydt, 

Jo Jr^ J \\x-Xo\\<R+M{To-t) 

limsupLs < XMi [ [ ||V(F-div/)(t,y,-)||L^dydt, 

£,ri,e^o Jo Jr'^ 

lim sup Lt = . 

£,f],e-*o 

Collating all the obtained results and using the equality ||V/i(x)|| = — pT+r/^'i 
[ [ \u{T,x) - u{T,y)\ -^/ii I j dxdy 

Jri^ J\\x-Xo\\<R+M{To-T) ^ V / 

< I I \u{0,x) -u{0,y)\^fii( ^ ^ ] dxdy 

Jr^ J\\x-Xo\\<B.+M(To-T) ^ V / 

-||V9«/||loo / / / \u{t,x) - u{t,y)\ 

Jo Jr'^ J \\x-Xo\\<R+M{To-t) 

1 / f\\x -y\ 



X;^lvTT'^i\^ ^ — j \\x-y\\dxdydt (4.21) 

+ {N\\Vduf W^oo + \\duF\\^^) r [ [ \u{t,x)-u{t,y)\ 

Jo Jr^^ J\\x-Xo\\<R+M(To-t) 

]^ 1 1 X — 1 1 \ 
X ^1 ( 1 dx dy dt 



+AMi r [ \\\/{F-dwm,y,-)\\^^dydt. 
Jo Jm^ 



If ||V9u/||loo = ||<9„F||loo = and under the present assumption that Uq & C 
using Proposition 14.31 (|4.8|) and (|4.2U|) . we directly obtain that 

TV(n(r)) < TVK) + ^ r / ||V(F-div/)(t,y,.)||Loodydt. (4.22) 

^1 Jo Jrn 

The same procedure at the end of this proof allows to extend (j4.22p to more general initial 
data, providing an estimate of TV (u{t)^ in the situation studied in [3J. 

Now, it remains to treat the case ||V5u/||loo ^ 0. A direct use of Gronwall type 
inequalities is apparently impossible, due to the term with V//. However, introduce the 
function 

T{T,X)= [ [ [ \u(t,x) — u{t,x — z)\ ^ IJ.1 f^^Y^] dxdzdt 

Jo Jr'^ J \\x-Xo\\<R+M{To~t) X V A 



SO that 



~T i i i x) — n(t, X — z)! ^^lyj^!^ ^ ||z|| dx dz dt . 

A Jo JM^ J\\x-Xo\\<R+M{To-t) ^ 

Denote C{T) = Mi [ [ \\V{F - divf){t,y, ■)\\ ^dy dt and integrate (ICTD on [0,r'] 
Jo Jm 

with respect to T for T' < Tq. It results 
^^(r',A) < ^/ / \u{0,x)-u{0,y)\fi{x-y)dxdy 

^ JRN J \\x-Xo\\<R+MTo 

+T' ||V9„/||loo dxHT', A) + y (2iV|| V5„/||l^ + \\duF\\^^) T{T' , A) 
+T' C{T'). 

Denote a = ( 2A^|| V5„/||loo + ||i9„F||loo — ^] {W^dufWi^^) ^, so that hmy/^o" = — oo. 



The previous inequahty reads, using (|4.1ip for Uq, 

dxHT'A)+a ^^\'^^ > -(MiTVK) + C(r')) ^ 



1 



5A(A"^(r',A)) > -A" (MiTV(uo) + C7(T')) 



|va„/|iLoo ■ 

Finally, if T' is such that a < —1, then we integrate in A on [A, +oo[ and we get 

i^(T', A) < ^— - (MiTV(no) + C{T')) \ . (4.23) 
A -a - 1 ' II Vd„/|lLoo 

Furthermore, by (14. 1|) and (14. 2p . there exists a constant -fC > such that for all z G 

-^;(||z||)<K^i(M^ . (4.24) 
Divide both sides in (|4.2ip by A, rewrite them using (|4.23p . (|4.24p . apply (|4.1ip and obtain 
v/ / \u(T,x) -u{T,y)\^Hi [^^^—-^) dxdy 

JRI^ J\\x-Xo\\<R+M{To-T) V ^ / 

'111 

2A 



< MiTV(no) + -^(^^ 2A) 2iv+2^||va,/||^^ + (2iV||V5,/||Loo + ||5„F||lo.) 



+Mi / / ||V(F-div/)(t,y,-)||Loodydt. 



An application of ()4.23p yields an estimate of the type 

V / 

' J B{xo,R+M{To-T)) 



A 



u{T, x) - u{T, x-z)\ ii{z) dxdz<C , (4.25) 



the positive constant C being independent from R and A. Applying Proposition 14.31 we 
obtain that u{t) £ BV(M^;M) for t G [0,2Ti[, where 

1 



2((l + 2iV)||Va„/||Loo + ||a„F||Loo) 



(4.26) 



The next step is to obtain a general estimate of the TV norm. The starting point 
is Km . RecaU the definitions (fCTjl of Mi and (IMj) of Ti. Moreover, by gS]), 



\\zr^l[i\\z\\)dz = -{N + l)M,. 

Divide both terms in (j4.2ip by A, apply (j4.9p on the first term in the right hand side, 
apply ()4.1ip on the second and third terms and obtain for all T S [0,Ti] with Ti < Tq 

TV(«(r)) < TVK) + ((2iV + l)||V5„/||Loo + ||5,F||Loo)^ / T\ {u{t)) dt 

^1 Jo 

+^ f f ||V(F-div/)(t,x,-)|L^dxdt. 
Jo iR^ 

An application of Gronwall Lemma shows that TV {u(t)) is bounded on [0, Ti]. Indeed, 

TV(n(t)) <e'^°*TVK) + ^ /"""e^-^^-*) [ \\V{F - dWf){t,x,-)\\^^dxdt (4.27) 

(-'I JO Jw^ 

forte [0,Ti], Mi,Ci as in Km . (O]) and = [(2iV + 1)|| Va„/|lLoo + ||a„F||Loo]Mi/C7i. 

We now relax the assumption on the regularity of Uq. Indeed, let Uo € BV(R^;]R) 
and choose a sequence of C^(M^;M) functions such that TV«) ^ TV(n o), as m 
Theorem 14.21 Then, by Theorem 12.31 the solutions li" to (II. ip with initial datum 
satisfy 

lim = u{t) in Li^„^ and TV (u{t)) < liminf TV (n"(t)) , 

n^+oo n— >+cxD 

where we used also the lower semicontinuity of the total variation. Note that (I4.27|) . as 
well as the relations above, holds for all t E [0,Ti], Ti being independent from the initial 
datum. Therefore, the bound (I4.27P holds for all BV initial data. 

Remark that the bound (j4.27p is additive in time, in the sense that applying it iteratively 
for times Ti and t yields ()4.27p for time Ti +t: 

TV (n(ri + 1)) 

< e^°'T\{u{Ti)) + ^ r^\^°^'~'^ f ||V(F-div/)(s,x,-)||Loodxds 
'^1 Jti Jrn 



Ml 



e'^-'^i TVK) + ^ [ \-^iT,-s) I ||v(F-div/)(s,x,-)||L.odxds 
t^l Jo JRN 



/'^'''*^"°^'"''"*"'^ / ||V(F-div/)(s,x,-)||Locdxds 

Ti+t 



Ci 

„Ko(Tl+t) 



TVK) + :^/"' e-o{T,+t-s) I ||v(F-div/)(s,x,-)||Loodxds, 
Jo Jr^ 



The bound ()4.27p can then be applied iteratively, thanks to the fact that Ti is independent 
from the initial datum. An iteration argument allows to prove (j2.5p for t G [0, To]. The 
final bound ()2.5p then follows by the arbitrariness of Tq, thanks to ()2.1|) . □ 



5 Proof of Theorem 12.6 



,-+), ^ € C 
pN 



The following proof relies on developing the techniques used in the proof of Theorem 12. 5i 

Proof of Theorem IMM Let $ G C^(M+ x M^;l+), ^' € C^^dR x ]R^;M+) and set 
(p{t, X, s, y) = ^{t, x)^{t - s,x -y) as in (|4.14p . 

By Definition EH we have V/ G M, V(t,x) G M+ x 



(5.1) 



(n - /) 9^9? + (/(s, y, u) - /(s, y, I)) ■ Vyf + (F(s, y, u) - div/(s, y, I)) (p 

xsign(u — l)dyds > 

and VA; G E, V(s, y) G M+ x 

- A;) dt(p + (^(t, X, f ) - g{t, x, k)) ■ V^ip + {G{t, x, v) - divg{t, x, k)) (p 

xsign(t; — k) dx dt > 0. 



(5.2) 



Choose k = u{s, y) in ()5.2p and integrate with respect to (s, y). Analogously, take / = v{t, x) 
in (|5.ip and integrate with respect to (t,x). By summing the obtained equations, we get, 
denoting u = u{s,y) and v = v{t,x): 



(u - v)'^dt^ + {g{t, X, u) - g{t, x, v)) ■ (V^>) 

+ {g{t, X, u) - g{t, X, v) - f{s, y, u) + f{s, y, v)) ■ (V5')$ 

+ {F{s, y, u) - G{t, X, v) + dwg{t, x, u) - div/(s, y, v)) 99 

xsign(n — i;) dx dt dy ds >0. 

Introduce a family of functions {^,?}^>o as in (|4.16p . Let M = \\dug\\i^oo(^^.^N^ and define 
X,i^ as in ()4.17p . for e,9,To,R > 0, Xq G M^, (see also Figured]). Remind that with 
these choices, equalities (j4.18p still hold. Note that here the definition of the test function 
(f is essentially the same as in the preceding proof; the only change is the definition of 
the constant M, that is now defined with reference to g. We also introduce as above 

the function B{t, x, u, v) = M\u — v\+ sign(ii — v) {g{t, x, u) — g{t, x, v)) ■ ^7 that is 

positive for all {t,x,u,v) G $7 x M^, and we have: 



< 



< 



{u — v)dt^ + {^g{t, X, u) — g{t, x, v)) • V<I> sign(u — v) dx dt dy ds 
/ / / / [\u - vlx'i' - B{t,x,u,v)xYQ]'^ dxdtdyds 

Jr+Jr^ Jr+ Jw^ 

/ / / / \u — v\ X 4^ dx dt dy ds . 

Jr4^Jrn Jr+Jr^ 



Thanks to the above estimate and (15.31). it results 



(n — f )xV^ 

+ {g{t,x,u) - g{t,x,v) - f{s,y,u) + f{s,y,v)) ■ (V^')$ 
+ {F{s, y, u) - G{t, X, v) + divy(t, x, u) - div/(s, y, v)) ip 
X sign(n — v) dx dt dy ds 



> 0, 



i.e. I + Jx + Jt + K + + Lt > 0, where 

I = \u - v\xip^ dxdtdyds , (5.4) 

Jr+ Jr^ Jr+ Jr^ 

J- = I I I I {f{t,y,v)- f{t,y,u) + f{t,x,u)-f{t,x,v)) -(V^)^ (5.5) 
Jr+ Jr^ Jr+ Jr^ 

X sign(n — v) dx dt dy ds , 

Jt = [ [ [ [ {f{s,y,v)-f{s,y,u) + f{t,y,u)-f{t,y,v))-{\/^)'^ 

Jr+ Jr'^ Jr+ Jr'^ 

X sign(n — v) dx dt dy ds , 

K = [ [ [ [ {{g-f){t,x,u)-{g-f)it,x,v))-{V^)<^ (5.6) 
Jr+ Jr^ Jr+ Jr^ 

X sign(n — v) dx dt dy ds , 

Lx = I / / / {F{t,y,u) - G{t,x,v) +d\Yg{t,x,u) - d\Yf{t,y,v)) if 
Jr+ Jr^ Jr+ Jr^ 

X sign(n — v) dx dt dy ds , 

Lt = / / / {F{s,y,u) - F{t,y,u) +divf{t,y,v) - dWf{s,y,v)) ip 

Jr+ jrJv Jr_^ Jrn 

X sign(n — v) dx dt dy ds . 



Now, we choose ^(t,x) = i^it) fi{x) as in (I4.19p . (liTT]) . (fO]) . Thanks to Lemma [52 
Lemma 15.31 and Lemma 15.41 we obtain 



hmsup/ < / \u{0,x) — v{0,x)\dx (5-8) 

e,'n,X-*0 J\\x-Xo\\<R+MTo+e 



\u{T, x) — v{T, x)\dx , 

\\x-Xo\\<R+MiTo-T) 

HmsupJ^ < N\\Vduf\\i^oc [ [ \v{t, x) - u{t, x)\ dx dt , (5.9) 

e,7?,A^0 Jo Jb{xo, 



,R+M(To~t)+b 



hm sup Lx < 

£,?7,A^0 Jo J B{xo,R+M{To-t)+e) 



dy dt 

Loo 



((F_G)-div(/-<7)) (t,y,-) 
+ ( iV||Va„/||Loo + ||5„F||loo + \\du{F - G)||loo) (5.10) 



X / / \v{t,x) — u{t,x)\dx6.t . 

Jo J B{xo,R+M{To-t)+e) 

Besides, we find that: 

\Jt\ < vWdtdufhoo [ [ [ [ \u{t,x)-u{s,y)\\\Vn^dxdtdyds, 
Jr+ Jr^ Jr+ jrJv 

\Lt\ < v^N{R + MTo)''{T + e){\\dtdwf\\j^^ + \\dtF\\^^), 
so that 

limsup I Jt[ = limsup |Lj| = . (5-11) 

r)-*0 ri^O 

In order to estimate K as given in (15. 6|) , we introduce a regularisation of the y dependent 
functions. In fact, let Pa{z) = (f ) and ^"^(y) = -JnCt y^j, where p € C^(M;M+) and 



jx.ivi. kjuiuuiuu, ivi. ivmivmi <y, ivi.u. jxusiiu 



a G C^(M^^;]R+) arc such that ||p||l1(]k.]k) ■ 
supp((t) C B{0, 1). Then, introduce 

P{w) = (g - f ){t,x,w) , Sa 
Tiiw) = Sc,{w - v) {Pi{w) - Pi{v)) , up 
T^w) = s\gn{w - v) {Pi{w) - Pi{v)) , 

so that we obtain 



^llLi(ffiJv.K) = 1 and supp(p) C ]-l,l[, 



sign *„ Pa , 



= / sign{w) {paiufs - V - w) Pi{up) - pa{u-v - w) Pi{u)) dy^(fdwdy 
- / sign(u;) [paiu/s - v - w) - pa{u - v - w)) Pi{v) dy^ipdwdy 

JR^ jR 

= [ [ r\ign{w)p'^{U -v-w){Pi{U)-Pi{v)) dy^ipdUdwdy 
Jr^ Jr Ju 

+ / / sign{w)pa{U -v-w)Pl{U)^y^ipdUdwdy. 

Jrn Jr Ju 

Now, we use the relation duSa{u) = ^p (^) to obtain 

(nM-n(u),a,,^) 

< / - sup ip(— — -] {PiiU) - Pi{v))] mm\2a,\u- up\] dy^ipdy 

Jm ot u&[{u,u0)] \ \ oi J y L ' 'J 

+ / / \Pl{U)\dy^^dUdy. 

JrN J a 

When a tends to 0, thanks to the Dominated Convergence Theorem, we obtain 

{T\up)-T\u),dy,<p) < [ \u-up\\\Pl\\^^dy,<pdy. 

Applying the Dominated Convergence Theorem again, we see that 

lim limiTJug), dy^ip) = {T{u), dy^ip) , 

lini lim(Ta(u^), V2;<^) = {T{u),Vyip). 

/J— »0 a— >0 

Consequently, it is sufficient to find a bound independent of a and P on Ka^i3, where 
Ka,f3 = - / / / Taiup) -Vyipdxdtdyds . 

Jr+ Jr^^ JR+ JRf^ 

Integrating by parts, we get 

Ka,i3 = BivyTa{ui3)(pdxdtdyds 



/ / / / duSa{ui3 -v)'Vui3- {{g - f){t,x,ui3) - {g - f)(t,x,v)) (pdxdtdyds 

Jr+ Jr^ Jr+ Jr^ 

/ / / / Sa{ui3-v){du{g- f)it,x,u/3)-Vup)(pdxdtdyds 

Jr+ JRf^ Jr+ Jffi^ 



+ 

= KI + K2. 

We now search a bound for each term of the sum above. 



For Ki, recall that duSa{u) = (^)- Hence, by Dominated Convergence Theorem, 
we get that ETi — > when a — > 0. Indeed, 



< 



2 f ur — V 

-P 

a \ a 

2 f Ur — V 

-p ' 



• {{g - f){t, X, Up) -{g- f){t, x, v)) ip 
ip\\Vup{s,y)\\ I \\duif - g){t,x,w)\\<lw 



a \ a ^ 

< 2||p||loo(ir;K) pup{s,y)\\ \\duU - g)\\^^^^.^M-^ip G x M^)2;m) . 

• Concerning K2, 

rT+e+ri r 

K2 < \\du{f-9)\\ 

i-T+e+rj 

< \\Mf - 9)Lo.^n;RN^ / TV(n^(t))dt. 

I' 

Finally, letting q, /3 — > and e, r/, A ^ 0, thanks to [H Proposition 3.7], we get 

limsupif < ||a„(/-5)||Loo /^TV(n(t))dt. (5.12) 

£,??,A-*0 JO 

Now, we collate the estimates obtained in ()5.8p . (15.91) . (I5.10|) . (15. lip and ()5.12p . Remark 
the order in which we pass to the various limits: first £,r],6 and, after, A ^ 0. 
Therefore, we get 



|n(r, x) — v{T, x) \ dx 



< 



B{xo,R+M{To-T)) 

|n(0, x) — v{0, x)\ dx 
B{xo,R+MTo) 

+ [2Ar||Va„/||Loo + \\duF\\^^ + \\du{F - G)||l. 



T 



+ 

or equivalently 
where 

A{T) = 



J B{xo,R+M{To-t)) 

\\duif - g)\\^,^ rTViuit))dt 



\v{t,x) — u{t,x) \ dxdt 








{{F-G)-dw{f-g)) {t,y,.) 

lO JB{xo,R+M{To-t)) 

A'{T) < ^'(0) +kA{T) + S{T) , 
\v{t, x) — u{t, x) \ dx dt , 



dy dt 



(5.13) 




'0 J B{xo,R+M{To-t)) 

K = 2N\\Vduf\\^o. + \\duF\\^^ + \\du{F-G)\\^^, 



(5.14) 



S{T) = \\dM-9)Lo. I TV(^.(t)) dt 

Jo 



+ 




J B{xo,R+M{To-t)) 



{{F -G)- div(/ - g)) (t, y, •) dy dt. (5.15) 

^ ' L°° 



The bound ([23]) on TV {u{t)) gives: 



< ^ a+ / ^ b{t)dt+ / c(t)di 

/^o JO Jo 



where Kq is defined in (j2.4|) and 

o = - 5)||LooTV(no) , 

6(i) = iVVF7v||a„(/-5)||Loo /" ||V(F-div/)(t,x,-)||L..dx, 



'B(a;o,_R+M(To-t)) 

since T < Tq. Consequently 



((F_G')-div(/-5)) 



(KoT 1 /"T Ko{T—i) 1 /"T \ 

° ~ a+ / —h{t)At+ / c(t)dt . (5.16) 
'to Jo Jo J 

By a Gronwall type argument, if = k, we get 



c{t) dt 



yielding 



/ \u{T,x) — v{T,x)\dx < e'""^ / \uo{x) — Vo{x)\ dx 

J\\x-Xo\\<R J\\x-Xo\\<R+MT 

+ Te^^TY{uo)\\du{f-g)\\^^ (5.17) 

+ NWn{ r {T-t)e<^~'^ f \\V{F-d\wf){t,x,-)\\^o^dxdt\\\du{f -g)\\^^ 

+ Te'^^^-*) / ((F_G)-div(/-5))(t,x,-) dxdt 

Jo J\\x-Xo\\<R+M(T-t) 

while, in the case Kq ^ k, we have 



a 



T pKo{T-t) _ „K.(T-t) 



Kq — H Jo 

Taking T = To, we finally obtain the result. 



T 



b{t)dt+ I e''(^-*)c(t)dt. 
'o 

□ 



Remark 5.1 Assuming that also {g,G) satisfies (H2), allows us to exchange the role of 
u and V in ^5A^. Let 



NWn [{2N + l)||Va„5|lL- + \\duG\\^^) 
\\duU -g)\\^^'I\{vo), 



m 



dM - 9) j^o^NWn / \7{G-divg){t,x,-) l^^dx, 
Jrn 

2iV||Va„(7|lL- + WduGWj^oo + \\du{F - g)||l^ , 



and repeating the same computations as above, we obtain 

A {T) < A{{)) + kA{T) + \ z a+ z b{t)dt+ c{t)dt 

\ '^o Jo '^o Jo ^ 

so that, finally. 



A'{T) < A'(0) + min(K, R)A{T) + max 



a + I 6(t)dt, 



z a+ / ; b{t)dt 

l^o Jo 



f^o Jo 
T 



+ 



[ c{t)dt. 
Jo 



We collect below some lemmas that were used in the previous proof. The first one 
reminds a part of the proof of [31 Theorem 2.1]. 



Lemma 5.2 Let I be defined as in |5.^[ j. Then, 

limsupl < / \u{0,x) — v{0,x)\dx 

e^O J\\x-Xo\\<R+MTo+8 

- \u{T,x) -v{T,x)\dx + 2 sup TV (u(t)) A 

J\\x-Xo\\<R+M{To-T) tG{0,T} 

+2 sup / \u{t,y) - u{s,y)\dy . 

te{0,T} J \\y-xJ\<R+X+M(To-t)+e 
se]t,t+vl 

Proof. By the triangle inequality I < h + h + h, with 

h = [ [ [ [ \u{t,x) - v{t,x)\x'{t)ij{t,x)^{t - s,x - y)dxdtdyds , 

\u{t,x) — u{t,y)\ |x'(0| "^{t — s,x — y) dxdtdyds , 

-^3 = /_ /_.,/_ I__\u{t,y)-uis,y)\\x'it)\ij{t,x)^{t-s,x-y)dxdtdyds 
Then, 

Ii = / / \u{t,x) -vit,x)\{Y^it) -Y^{t-T))'iP{t,x)dxdt 



< \u{t,x) - v{t,x)\Y^{t)dxdt 

Jr+ J \\x-Xo\\<R+M {To-t)+9 

- / \u{t,x) -v{t,x)\Y^{t-T)dxdt 

Jr+ J \\x-Xo\\<R+M{To~t) 
and by the L-*- right continuity of u and « in time, thanks to Theorem 12.31 

limsup/i < / \u{0,x) — v{0,x)\dx 

e^O J\\x-Xn\\<R+MTo+e 

\u{T, x) — v(T, x)\dx . 

x-Xo\\<R+M{To-T) 



For I2 and I3, we have 
I2 < [ [ [ |u(t,x)-n(t,y)|(y;(t) + y;(t-r))/xdxdydt, 

Jr+Jr^ J\\x-Xo\\<R+M{To-t)+e 



h < \u{t,y)-u{s,y)\{Y^{t)+Yl{t-T))vdydsdt. 

JR+JR+J\\y-Xo\\<R+\+M{To-t)+e 

As e — i- 0, we use on the one hand the L-*- right continuity in time of u, thanks to The- 
orem [231 and on the other hand that u{t) E BV(M^;M), thanks to Theorem 12.51 In 
particular, we can use (|4.1ip to obtain 

hmsup/2 < 2^ sup / \u{t,x) — u{t,x h)\Ax 

e^O l|/i||<A J \\x-Xo\\<R+M(Ta-t)+e 

< 2 sup / \u{t,x) — u[t,x + h)\dx 

\\h\\<\ J\\x-Xo\\<R+M(To-t) + e 
t6{0,T} 

< 2 sup TV (n(t)) A , 

te{o,T} 

hmsup/3 < sup / \u{t,y) -u{s,y)\dy 

rp s&]t,t+rj[J \\y-Xo\\<R.+\+M{To-t)+e 



t=0 



< 2 sup / \u{t,y) - u{s,y)\dy . 



te{o,T} J \\y-x„\\<R+X+M{To~t)+e 

se]t,t + r)[ 



□ 



Lemma 5.3 Let Jx be defined as in ( 15.5)) . Then, 



hmsupJa: < N\\Vduf\\i^oo I I \v{t,x) — u{t,x)\dxdt 



I B{xo,R+M{To~t)^ 

+NT\\Vduf\\^oo sup T\{u{t))\ 



+NT\\Vduf\\^^ sup / 

*6[0,T1 J\\y-Xo 



u{t,y) - u{s,y)\ dy. 



\\y-Xo\\<R+X+A4{Tc-t)+e 

se]t,t+r]l 

Proof. By assumptions (HI), / G C^(0;M^) and therefore 

\\f{t,y,v) - f{t,y,u) + f{t,x,u) - f{t,x,v)\\ 

/'v(t,x) fl 

/ / \/duf {t,x{l — r) + ry,w) ■ {y — x) dr dw 

'u{s,y) Jo 

< I|V9„/||loo \\x - y\\ \v{t,x) - u{s,y)\ . 



Then, 



Jx < ||V5„/||loo / / / / \vit,x) -u{s,y)\\\x-y\\\\Vfi\\iyxi^dxdtdyds. 
Jr+ jr^ Jr+ jr^ 




IR+ JR^ JR+ 



Similarly to the proof of Lemma \5.2\ we apply the triangle inequality and obtain < 
Ji + J2 + J3 where 

Ji = IIv5„/||l^ 

J2 = ||Vd„/i|Loc 

J3 = I|V5„/||loo 
For Ji, we have, thanks to (|4.5|) 



Ji < N\\vduf w^oo r^' [ 

Jo J B 



\v{t, x) 


— u{t, x) 


\\x 


-y\\ l|V/u||i. 


X V' da; dt dy ds 


|n(t, x) 


- u{t,y) 


\\x 


-y\\ liv^||i^ 


Xtpdxdt dy ds 


\u{t,y) 


- u{s,y) 


\\x 


-y\\ llV/i|li^ 


Xipdxdt dy ds 


{xo,R+M(To~t)+e) 


v{t 


x) — u{t, x) 


dx dt. 



For J2, we have 



rT+e r 

J2 < iV||Va„/||Loo / sup / \u{t,x) - u{t,x + h)\dxdt 

Jo \\h\\<\ J \\x-Xo\\<R+M {To-t)+e 

< N\\Vduf\\^o.{T + e) sup TY{u{t))X, 

TG[0,T+e] 



and for J3 



J3 < iv||va„/||L^ / ' sup / 

Jo se]t,t+-q[J\\x-Xo 

< iv||va„/|iLoc(r + £) sup [ 

te[0,T+el J\\y-Xo 



\\<R+\+M{To-t)+e 
y-Xo\\<R+X+M(To -t)+e 



u{t,y) - u{s,y)\ dydt 
u{t,y) - u{s,y)\ dy. 



se]t,t+ril 

In particular, letting A, r],e,6 ^ 0, we prove that ^2, <^3 ^ and 

rT 

\v{t,x) — u(t,x) \ dxdt 



limsup Ji < iV||Va„/||Loo [ [ 



R+M{To-t)) 



completing the proof. 



□ 



Lemma 5.4 Let he defined as in (5.1) and Mi as in ^.20 ). Then 



lim sup Lx < T 

e^O Jo J\\x-Xo\\<R+M(To-t)+6 



((i7_G)-div(/-5)) (t,x,-) 



dx dt 



+ [N\\Vduf\\^^ + \\duF\\^^ + \\du{F - 

' / \v{t,x) — u{t,x)\dxdt 

J B{xo,R+M{To-t)+e) 

+T sup TV(u(t))A 

Te[o,T] 



+T sup / \u{t,y) - u{s,y)\dy 

t6[0,Tl J\\y-x„\\<R+\+M{To-t)+e 
se]t, t+77[ 

+AMi r [ ||V(F-div/)(t,x,-)||Loodxdt. 
Jo Jrn 



jx.ivi. kjuiuuiuu, ivi. ivmivmi <y, ivi.u. jxusiiu 



Proof. Let 

Li = / / / {{F -G) -divif - g)) {t,x,u)(psiga{u-v)dxdtdyds , 

Jr^ .7R+ Jr^ 

L2 = / / / {{F -G){t,x,v) - {F -G){t,x,u)) ipsign{u-v)dxdtdyds, 

Jr+ jr^ Jr+ jr^ 

^3 = / / / / {Fit,y,u)-Fit,y,v)+diyf{t,x,u)-diyf{t,x,v))^ 
Jm.+ j]r^ Jm.+ j]r^ 

X sign('u — v) dx dt dy ds , 

L4 = / / / / {{F -dwf){t,y,v) - {F -dwf){t,x,v)) (psiga{u-v)dxdtdyds, 
Jr+ Jrn Jr+ Jw.^ 



so that Lx = Li + L2 + + L4. Clearly, 

Jo J\\x-Xo\\<R+M(To-t)+e 

For L2 and L3, we have 



((G-F)-div(/-5)) 



dx dt . 



< ||5„(F-G)||j^^ / [ [ f \uis,y)-vit,x)\<pdxdtdyds, 

Jr, Jrn Jr, Jrn 




sign(u — I y (^dudiy f {t,x,w) + duF{t,y,w)) dwj (pdxdtdyds 



< (iV||Va„/||Loo + ||5„F||loo) / Iff \v{t,x)-u{s,y)\ipdxdtdyds. 

Jr+ Vffi^v 7]Rjv 

Proceeding as for J^, we find the following bound for //// x) — u{s,y)\ip in L2, L3. 
i:2 + i^3 < (Ar||V5„/||Loc + ||a„F||Loo + ||5„(F-G)||j^^) 



T+ 



\v{t, x) - u{t, x) I dxdt + (T + e) sup TV {u{t)) X 

^B(j;o,-R+M(ro-t)+6>) T6[0,T+e] 



+(r + e) sup / 

teIO,T+e] dilu- 



te IO,T+e] J||j,_a;„||<fl+A+M(To-t)+6» 

sS]t,t+r)[ 



\u{t,y) - u{s,y)\ dy 



For L4 we have 

U = [ [ [ [ [ V{F-diYf){t,rx + {l-r)y,v)-{y-x)dr 

Jr+ Jrn Jr+ Jr^ [jo 

X sign(M — i;) dx dt dy ds 
< AMi / ||V(F-div/)(t,x,-)||Locdxdt. 

To complete the proof, it is sufficient to note that = Li + L2 + + L4. 



□ 



jx.ivi. kjuiuuiuu, ivi. ivmivmi <y, ivi.u. jxusiiu 
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